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ABSTRACT 

We present the discovery of the Kepler-19 planetary system, which we first identified 
from a 9.3-day periodic transit signal in the Kepler photometry. From high-resolution 
spectroscopy of the star, we find a stellar effective temperature T e ff=5541 ± 60 K, a 
metallicity [Fe/H]=— 0.13 ± 0.06, and a surface gravity log(g)=4.59 ± 0.10. We combine 
the estimate of T e g and [Fe/H] with an estimate of the stellar density derived from 
the photometric light curve to deduce a stellar mass of M* = 0.936 ± 0.040 M and 
a stellar radius of R± = 0.850 ± 0.018 Rq (these errors do not include uncertainties 
in the stellar models). We rule out the possibility that the transits result from an 
astrophysical false positive by first identifying the subset of stellar blends that reproduce 
the precise shape of the light curve. Using the additional constraints from the measured 
color of the system, the absence of a secondary source in the high-resolution spectrum, 
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and the absence of a secondary source in the adaptive optics imaging, we conclude 
that the planetary scenario is more than three orders of magnitude more likely than 
a blend. The blend scenario is independently disfavored by the achromaticity of the 
transit: we measure a transit depth with Spitzer at 4.5 fim of 547^}^ ppm, consistent 
with the depth measured in the Kepler optical bandpass of 567±6 ppm (corrected for 
stellar limb-darkening). We determine a physical radius of the planet Kepler-19b of 
R p = 2.209 ± 0.048 i?®; the uncertainty is dominated by uncertainty in the stellar 
parameters. From radial-velocity observations of the star, we find an upper limit on 
the planet mass of 20.3 M®, corresponding to a maximum density of 10.4 g cm~ 3 . 
We report a significant sinusoidal deviation of the transit times from a predicted linear 
ephemeris, which we conclude is due to an additional perturbing body in the system. We 
cannot uniquely determine the orbital parameters of the perturber, as various dynamical 
mechanisms match the amplitude, period, and shape of the transit timing signal and 
satisfy the host star's radial velocity limits. However, the perturber in these mechanisms 
has period < 160 days and mass < 6Mj up , confirming its planetary nature as Kepler- 
19c. We place limits on the presence of transits of Kepler-19c in the available Kepler 
data. 

Subject headings: eclipses — stars: planetary systems — stars: individual (Kepler-19, 
KOI-84, KIC 2571238) 



1. Introduction 



With the recent discoveries of the first transiting ex oplanets int ermediate in size betwee n 



Earth and Neptune, namely CoRoT-7b (ILeger et al.l 12009). G J 1214b (ICharbonneau et al.l 120091 ) 



Kepler-9d ([Torres et al.ll201lh . Kepler-lO bc (IBatalha et al 



201ll : lFressin et al 



2011 



, Kepler-llbcdf 

(jLissauer et al.ll201ll ). and 55 Cancri b (jWinn et al.ll201ll ). astronomers have begu n in earnest to 



probe this radius regime of exoplanets, for which no Solar System analog exists. iBorucki et al 
(|201ll ) presents a catalog of 1235 transiting planetary candidates, of which nearly 300 have a 
radius estimate in the range 1.25< R v <2.0 R ffi . Whil e most of these candidates have not yet 
been confirmed as authentic planets, iMorton Johnson! (|201ll ) have shown that the rate of false 
positives is expected to be low for the Kepler-identified sample. T he composition of such pl anets 
may be widely variable, as exemplified by the case of GJ 1214b (jCharbonneau et al.1 12009). for 
which the measured radi us and mass were consistent with both a hydrogen envelope or a pure 
CO"2 or H2O atmosphere (jMiller-Ricci &: Fortnevil2010l ). It was only with follow-up studies of the 
atmosphere in transmission that it became possible to distinguish among the various possibilities 



for the composition of GJ 1214b (|Bean et al. 



2010 



Croll et al 



2011 



Desert et al. 



ous p o 
201 j ). 



The limiting precision of the current state-of-the-art radial velocity observations (meters per 
second) presents a challenge for the dynamical confirmation of these small planets. In the case of the 
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1.42 i?0 transiting planet Kepler- 10 b, iBatalha et al.l (|201ll ) gathered 40 high-resolution spectra at 
the Keck telescope ( Vogt et al.lll994l ) to determine a mass of 4. 56^29 M ffl . In the absence of radial 
velocity confirmation, however, it is still possible to make a statistical argument for the planetary 
nature of the candidate, if the combined likelihood of all false positive scenarios (namely, blends of 
stars containing an eclipsing member) is sufficiently smaller than the planet scenario. This process 
of "validation" for Kepier-identified plan etary candidates ha s already be en applied to three p lanets 

and 



Kepler-lOc (jFressin et al 



2011 



in the 1.5-3 radius range: K epler-9d (jTorres et al.ll201ll ). Kepler- 1 If (jLissauer et al 



2011 



)• 



Transiting planets are also of interest as they present an opportunity to identify yet more 
planets in the system by the method of transit timing variations (TTVs). Since this method 
was proposed (jHolman fc Murravl 120051 : lAgol et al.ll2005h. the sear ch for planets by TTVs has 
been a major activity in exoplanet research. ISteffen fc Agoll (|2005l ) applied a lack of significant 
TTV variations in the TrES-1 system to deduce constraints on th e exist ence of additional, non- 



transiting planets. S ubsequent w orks, such as iMiller-Ricci et al.1 (|2008l ) fo r HP 189733 , using 
MOST observations, beanl (|2009h for C0R0T-I, using the C0R0T satellit e. iGibson et al.l (12009 ) 



for TrES-3, using observations gathered at the Liverpool Telescope, and iBallard et al.l (|2010bl ) 
for GJ 436, using EPOXI observations, have also used transit times to rule out companions, 
specifically companions in resonances, for which the TTV method is particularly sensitive to low- 
mass planets. The Kepler team has pre sented two case s of transit timing va riations in exoplanetary 



systems: Kepler-9 (jHolman et al.ll2010l ) and Kepler-11 (jLissauer et al.ll201ll ). In both of these cases 



the additional planet (or planets) responsible for the transit timing signal also transit, which enabled 
mutual constraints on the masses of the planets as predicted bv lHolman &: Murravl (|2005l ). Over 
the past year, the Kepler team has also presented instances of single-transiting candidate systems 
showing trans it timing variatio ns, but has not confirmed the planetary nature of the candidates 



or perturbers (jFord et al.ll201ll ). Meanwhile, using ground-based observations, several groups have 
described their transit times as being inconsistent with a linear ephemeris, though two such claims 
have been revisi ted b y groups who could not confirm the result. In the case of HAT-P-13, while 
Pal et al.l (|201ll ) and iNascimbeni et al.l (120111 ) found evidence for a companion from the transit 
times of HAT-P-13b, I Fulton et al.l (|201ll ) demonstrated that the times were consist ent with a 



linear ephemeris (with the exception of a single transit). In the case of OGLE-lllb, iDiaz et al 



(200 8]) claimed that the transit times were inconsistent with a linear ephemeris, but an analysis by 
Adams et al.l (j2010l ) with additional transit observations found no e yidence for TTVs o r dur ation 
variations a n d poi nted to systematic errors in previous photometry. IMacieiewski et al.l (j2010l ) and 
Fukui et al.l (|201ll ) presented evidence for transit timing variations in the WASP-3 and WASP-5 
system s, respectively, bu t cautioned that additional transits are necessary to con firm or refute the 



signal (jFukui et al 



2011 



expressed caution about unknown systematic effects). IMacieiewski et al 



(j201ll ) presented evidence for TTVs of WASP-lOb, and they reported a two-planet orbital solution 
that fit the TTVs and radial velocities better than alternative orbital models they found, which 
was not achieved in prior work. 



-4- 



In this paper, we present the discovery of two planets orbiting Kepler-19. The star, which has 
right ascension and declination 19h21m40.99s and +37d51m06.5s, Kepler magnitude Kp=11.90, 
and Kepler Input Catalog number 2571238, was identifi ed to host a p l aneta ry candidate in the 
catalog of 1235 Kepler identified candidates published bv lBorucki et al.l (|201ll ). In that work, the 
star was identified by the Kepler Object of Interest (KOI) designation KOI-84. The first planet, 
identified by its transits, has a period of 9.3 days and a radius of 2.2 Rq, as we discuss below. 
We validate the planetary nature of the transit signal by a blend analysis. The second planet, 
Kepler-19c, was identified by transit timing variations. We see no evidence for transits of Kepler- 
19c in the available Kepler data. This detection differs from the ones using ground-based data, 
summarized in the previous paragraph, in several ways. In the case of Kepler-19, the transiting 
object has a radius of only 2.2 R®, whereas other claims are for perturbations to the transit times 
of hot Jupiters. Additionally, as we show below, we have well-sampled the TTV signal, since we 
have measured the transit times of Kepler-19b at a cadence 30 times shorter than the TTV period. 
Finally, the TTV signal reported here is a much higher signal-to-noise detection. In Section 2, we 
present the Kepler time series from which we detected the system. In Section 3, we present our 
characterization of the Kepler-19 system from the photometry. And in Section 4, we summarize our 
follow-up observations of the star. In Section 5, we present the validation of Kepler-19b as a planet. 
In Section 6, we discuss our constraints on the nature of the perturbing planet Kepler-19c from 
transit timing variations, as well as our inferred constraints on the composition of the transiting 
planet Kepler- 19b. 



2. Kepler Observations 



The Kepler spacecraft, launched on 7 Mar ch 2009, will photomet r ically mon itor 170,000 stars 



for 3.5 years for evidence of transiting planets. IVan Cleve C aldwell (2009) and Argabright et al 



(2008) give an overview of the Kepler instrument, and 



Caldwell et al 



(|2010h and 



Jenkins et al 



2010bl ) provide a summary of its performance since launch. The Kepler observations of Kepler-19 



that we present in this work were gathered from 5 May 2009 to 5 March 2011. This range spans 
Kepler "Quarters" 0-8; Kepler operations are divided into four quarters each year. At the end 
of each quarter, Kepler rotates the spacecraft by 90° to maintain illumination of the solar panels 
(though the Quarter pre-science commissioning data were gathered in the same configuration as 
Quarter 1). For Quarters 3-8, the observations of Kepler-19 were gathered continuously with an 
exposure tim e of 58.8 seconds, corre spon ding to the "sho r t cade nce" (SC) mode of the instrument, 
described by iJenkins et al.l (|2010bl ) and iGilliland et al.l (|2010l ). while data from Quarters 0, 1, 
and 2 were gathered in long-cadence mode (characterized by an exposure time of 29.5 minutes). 
The data contain gaps of approximately 3 days between quarters for scheduled downlinks. We 
used the raw light cu rves generated by the Kepler aperture photometry (PA) pipeline, described in 
Twicken et al.1 (|2010l ). to which we add two additional steps. First, we remove the effects of baseline 
drift by individually normalizing each transit. We fit a line with time to the flux immediately before 
and after transit (specifically, from 9 hours to 20 minutes before first contact, equal to 2.5 transit 
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durations, and an equivalent time after fourth contact). We divide this line from the observations 
spanning five transit durations and centered on the predicted transit time. For the observations 
gathered outside of transit, we apply a median filter, with width equal to one day, in order to 
remove baseline drift over timescales of days. We observed slight flux offsets between observations 
that occurred after gaps of larger than 1 hour (either for data download, quarterly rolls, or safe 
modes). The Kepler time series of Kepler-19 are shown in Figure 1. 
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Fig. 1. — Kepler-19 relative photometry from Quarters 0-8, as a function of barycentric Julian 
Day. Data from quarters 0, 1, and 2 were gathered in long-cadence observing mode, while data 
from quarters 3-8 were gathered in the short-cadence mode. Top panel: The raw flux depicted in 
red (with bin size of 1/3 days). We have multiplied the short-cadence observations by a factor of 
30 to account for the exposure time ratio between modes, so that all observations appear on the 
same scale. The gaps between quarters are depicted by dashed lines. To remove the flux offsets 
between quarterly rolls, we compare the mean brightness during the first two hours of each quarter 
to the mean of final two hours of the previous quarter, and divide this ratio value from the flux 
over the entire quarter. This corrected flux is shown in black. Bottom panel: the detrended Kepler 
light curve, after applying a median filter with width equal to one day and normalizing individual 
transits as described above. This light curve is depicted with a binsize of 1 hour, so that individual 
transit events are apparent. Gaps of longer than one hour in the observations are associated with 
instances of correlated noise in the corrected light curve. 
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3. Analysis 

3.1. Derivation of Planetary Parameters from the Kepler Light Curve 

The traditional procedure for fitting the phased transit light curve relies on the assumption of 
a linear ephemeris. In this case, that assumption does not hold: the transit times deviate from a 
linear ephemeris in a nearly sinusoidal manner, with a period of about 316 days and an amplitude 
of about 5 minutes (as compared to a transit duration of 3.5 hours). We incorporated the transit 
timing devia tions into the l i ght cu rve parameter fit in an iterative sense, similar to the procedure 



described bv lLissauer et al.1 (|201ll ) for the transit times of the Kepler-11 planets: we first estimate 
the light curve parameters with an assumption of a linear ephemeris. We then fit for the epoch 
of each individual transit using these parameters in a manner that we describe below (with the 
exception of the epoch, which we allow to float for each transit). Finally, we shift the transits 
by their measured timing deviation, refold the light curve, and repeat the fit for the parameters. 
We repeat this process until all parameters converged, and found that it converged after only two 
iterations. We note that we employed only the subset of short cadence observations for fitting the 
light curve parameters themselves. The addition of the long cadence observations from Quarters 
0-2 to the analysis would have contributed only a very modest improvement to our knowledge of 
the parameters. 

We estimated the uncertainty in the parameters using the Markov Chain Monte Carlo (MCMC) 



metho d as follows, using a Metropolis-Hastings algorithm with G ibbs sampling iTegmark et al 



(|2004l ). We employ model light curves generated with the routines in lMandel fe Agoll (J2002J), which 
depend upon the period P, the epoch T c , the planet-to-star radius ratio R p /R+, the ratio of the 
semi- major axis to the stellar radius a/i?*, the inclination of the orbit i, and two quadratic limb- 
darkening coefficients u\ and U2- We assume an eccentricity of zero for the orbit, which we discuss 
in further detail in Section 3.2. We choose randomly one parameter, perturb it, and evaluate the 
y 2 of the solution. If the y 2 is lower, we accept the new parameter value. If the x 2 is higher, 
we evaluate the probability of accepting the jump as p = e"^* 2 / 2 . We adjust the width of the 
distribution from which we randomly draw the jump sizes in each parameter until 20-25% of jumps 
are executed in each of the parameters. We created five chains, each of length 10 6 points, where 
each of the chains is begun from a different set of starting parameters (each parameter is assigned 
a starting position that is +3<r or -3a from the best-fit values). We discard the first 20% of jumps 
from each chain to remove the transient dependence of the chain on the starting parameters. 

After the first iteration of the MCMC procedure, we locate the best solution in the y 2 sense, 
and use these parameters to fit the individual times of transit for observations gathered in both 
short and long cadence observing mode. We fix all parameters with the exception of the center of 
transit time, which we allow to vary over a range of 30 minutes for each transit, centered on the 
predicted transit time from a linear ephemeris, in intervals of 6 seconds. In this case, we account for 
the finite integration time by taking a numerical average of an oversampled model (evaluated at 6 
second intervals) over a period corresponding to the exposure time (58.8 seconds for short-cadence 



-7- 



observations, and 29.5 minutes for long-cadence observations). We determine the center of transit 
time from the epoch associated with the minimum x 2 value, and the error from the range over 
which Ax 2 < 1 when compared to the minimum value. In practice, this results in asymmetric error 
bars for many of the individual transits. We also conduct an analysis to determine the contribution 
of quarter-by-quarter correlated noise to the transit time measurements, whereby we inject transits 
with the same light curve parameters as we derive for Kepler- 19b into the Kepler light curve, and 
then fit for the transit times in an identical fashion to the authentic transits. Based upon a typical 
deviation of these times from the injected time (as compared to their error bars), we inflate the 
error bars for the short-cadence observations by factors of 1.32, 1.15, 1.45, 1.40, 1.18, and 1.20 for 
Quarters 3-8, respectively. These values are consistent with the larger scatter of the transit times 
in the latter quarters, as discussed in further detail in Section 6.1. The transits gathered in long 
cadence from Quarters 0-2 show errors that are consistent with Gaussian noise, and so we do not 
apply a scaling factor to these error bars. After two iterations of the steps described above, we 
found that the measured individuals transit times varied by less than 10 seconds between the two 
iterations. We fixed the average period to the one determined by the final fit to the transit times 
for the final MCMC analysis. 

In Figure [21 we show the MCMC correlations between all free parameters in the model fit, as 
well as the histograms corresponding to each parameter. In Figure [31 we show the Kepler transit 
light curve for Kepler-19b (which is phased to the best period after shifting the transit times by 
the values given in Table ED, with the best-fit transit light curve overplotted. We report the best-fit 
parameters and uncertainties in Table [TJ The range of acceptable solutions for each of the light 
curve parameters (R p /R+, a/R+, and i) is determined as follows. We report the "best" solution 
from the set of parameters that minimize the x 2 ■ The error bars are then given by the highest and 
lowest values that are within the 68% of points closest in \ 2 to the best value. We additionally 
ca lculate the transit duration, imp act parameter, a nd ingress duration, using the formulae given 
in ISeager &: Mallen-Ornelasl (|2003l ) and IWinnl (|2010l ) to create the distribution in those quantities 



from the parameters in the MCMC chains. In some cases, the error bar is asymmetric (for a/R* and 
i, which we expect from their asymmetrical MCMC distributions). We report the transit times, 
deviation from a linear ephemeris, and errors in Table [2J Figure [H shows the individual timing 
deviations from the best linear ephemeris. In Figure El we show a binned subset of "late" and 
"early" transits, comprised of five transits each (corresponding to numbers 26-30 and numbers 39- 
43 listed in Table [2]), over which we plot a model generated with a linear ephemeris. The deviation 
of the transit times from the predicted T c of five minutes, which is equal to approximately one 
ingress time, is apparent. 

We also performed the MCMC analysis with the two quadratic limb-darkening coefficients 
(LDCs), ui and 112, fixed to theoretical values. We used the values for the effective stellar temper- 
ature T e fj, metallicity [Fe/H], and surface gravity log(g) derived from our Spectroscopy Made Easy 
(SME) analysis of the Keck HIRES spectra (see Section 4.2): T cfr =5541 ±60 K, [Fe/H]=-0.13±0.06 
dex, and log(g) of 4.59 ±0.1 dex. In this case, the closest stellar model from the tables of theoret- 
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Fig. 2. — Markov Chain Monte Carlo probability distributions for light curve parameters of Kepler- 
19b. The dark grey area encloses 68% of the values in the chain, while the light grey area encloses 
95% of the values. We assign the range of values corresponding to la confidence from the area 
enclosing 68% of the values nearest to the parameters associated with the minimum x 2 (as described 
in the text). 
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Fig. 3. — Kepler transit light curve for the short-cadence observations of Kepler-19, centered on 
time of transit, with transit timing variations removed, and binned by a factor of 40 (binsize of 4 
minutes). Overplotted in red is the best transit model light curve, with parameters given in Table 
[TJ The bottom panels shows the residuals of the light curve, after the model is subtracted. 
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Fig. 4. — Kepler transit times for Kepler-19b from Quarters 0-8, as compared to the best linear 
ephemeris model. The linear ephemeris we use to generate these O-C values is given in Table [TJ 
and the individual transit times are given in Table [2j The demarcation between long cadence and 
short cadence observations is shown with a dotted line. 
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Fig. 5. — A binned subset of five "late" transits of Kepler-19b (in red, comprised of individual 
transit numbers 26-30 listed in Table [2]) and five "early" transits (in blue, comprised of transit 
numbers 39-43), with best linear transit model overplotted. At left, we show the binned transits 
centered on time of ingress, and at right, centered on time of egress. 




ical limb-darkening coefficients generated for the Kepler bandpass by iPrsal (|2010l ) corresponded to 
a model with T e fj=5500 K, log(g)=4.5, and [Fe/H]=0.0: these coefficients are «i=0.5 and U2=0.18. 
Our results for the planetary parameters were consistent with the values we obtained while allowing 
the limb-darkening coefficients to float, but the error bars were slightly smaller with fixed LDCs. In 
particular, the LDC derived from the light curve lie within la of the theoretical values for u\ and 
U2, and the largest deviation in the derived parameters is 1.4a for R p /R*. The difference between 
the best x 2 values (between fixing the LDCs or allowing them to float) is approximately equal to 
five, which is roughly consistent with the addition of two degrees of freedom. 



3.2. Physical Parameters 



We based our p rocedure for c onstr aining the mass, radius, and age of the host star on the 
method described bv lTorres et al.l (|2008l ). Using the metallicity determined from SME (described 
in Sectio n 4.2), we crea ted a set of stellar evol ution models from the Yonsei-Yale (Y 2 ) isochrone 
series bv lYi et al.l (|200ll ). with corrections from iDemarque et al.l (|2004j). We employed the interpo- 
lation software that accompanied that work, which accepts as inputs the age of the star, the iron 
abundance, and the abundance of a-elements relative to solar (for which we assume the solar value), 
and outputs a grid of stellar isochrones corresponding to a range of masses. We evaluated a set of 
isochrones over an age range of 0.1 to 14 Gyr (at intervals of 0.1 Gyr) and in [Fe/H] in 60 equally 
spaced increments from — 3a to +3a above and below the best-fit value of [Fe/H]=— 0.13 ± 0.06. 
We then performed a spline interpolation of each output table at a resolution of 0.005 Mq in ef- 
fective temperature T e g, the log of the surface gravity log(g), and the stellar luminosity L+. We 
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evaluate the physical radius corresponding to each stellar model via log(g) and the mass of the 
star (g = GM+/R 2 ), though it is also possible to convert to physical radius using the model stellar 
luminosity and effective temperature (assuming L+ = 47ri?*o"T 4 ); in practice these conversions give 
identical results. 



Re arranging Kepl e r's ve rsion of Newton's th i rd law in the manner employed bv lSeager &: Mallen-Ornelas 



(|2003l ). ISozzetti et al.l (120071 ) and iTorres et al.l (|2008l ). we convert the period (derived from pho- 
tometry), and the radius and mass of the host star (from isochrones) to a ratio of the semi-major 
axis to the radius of the host star, a/i?*: 



a 

R* 



Q \V3 p2/3 



4tt 2 



R* 



(M* + M % 



a/3 



(i) 



where we assume that M p is negligible when compared to the mass of the host star, and that the 
orbit is circular. Using the MCMC sequence of a/i?* and generating a series of Gaussian random 
realizations of [Fe/H] and T e g using the values and error bars derived from spectroscopy, we locate 
the best isochrone fit at each realization using the \ 2 goodness-of-fit 



X 



(J. 



a/R* 



+ 



cff 



V °T cff 



+ 



A [Fe/H] 



[Fe/H] 



a 



(2) 



Using the output of the MCMC chain of a/R± ensures that correlations between parameters, which 
are preserved in the chain, are properly incorporated into our estimate of the stellar parameters. 

We then assign a weight to the likelihood of each stellar model i n the chain by a pplying a prior 
for the initial mass function (IMF) that assumes a Salpeter index (|Salpeterlll955l ). The number 
of stars of e a ch ma ss and age, per 1000 stars, is generated by the interpolation software provided 
by lYi et al.l (|200ll ) for several IMF assumptions, including the Salpeter IMF. We designate the 
weight assigned to each stellar model in the chain by normalizing to the highest IMF value within 
the sample: in practice, the weights vary from 0.2 to 1 (from the least to most likely). We then 
incorporate this likelihood by discarding members of the chain according to their weight, where the 
weight is equal to the likelihood of remaining in the chain. About 50% of the original chain remains 
intact after this stage. The value for each stellar parameter is then assigned from the median of this 
weighted distribution, with the formal error bars assigned from the nearest 68% of values above 
and below the median. We find M*= 0.936 ± 0.040 M , R i! = 0.850 ± 0.018 R Q , and an age = 
1.9 ± 1.7 Gyr. These uncertainties exclude possible systematic uncertainties in the stellar models. 
Using the modified MCMC chain in both stellar radius and R p /R± to determine the physical radius 
of the planet, we find i? p =2.209 ± 0.048 Rq. The IMF prior changes the final answer by less than 
la for all physical parameters, and by O.ler in the case of the planetary radius. 

We note that we also recorded the log(g) of the best stellar model for each realization of 
T e ff, [Fe/H], and a/R±, as described above. From this analysis, we find a log(g) of 4.54±0.02, 
which is consistent with the value inferred from spectroscopy of 4.59±0.10. We conclude that 
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the assumption of zero orbital eccentricity for the transiting planet is consistent wit h the value of 



log[g) measured spectroscopically. However, using the analytic formulae presented in lCarter et al 



(200 8|), the derived log (g) would vary by only 0.04 dex if the eccentricity were as high as 0.15 (which 
found was typical for the sample of Kepler Objects of Interest), which is well 



Moorhead et al. 2011 



below our measured uncertainty on log(g). 
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Table 1. Star and Planet Parameters for Kepler-19 



Parameter 


Value 


Kepler-19 [star] 


Right ascension 


19h21m40.99s 


Declination 


+37d51m06.5s 


[Fe/H] 


-0.13±0.06 


log(g) [cgs] 


4.59±0.10 


T cS [K] 


5541±60 


v sin i [km s _1 ] 


< 2 


M* [M ] 


0.936±0.040 


-R* [i? Q ] 


0.850±0.018 


Age [Gyr] 


1.9±1.7 



Kepler-19b 



Period [days] 
T c [BJD-2450000] 
Rp/R* 
a/R* 
i [degrees] 

Ul 
U2 

Impact parameter, b 
Transit duration, r [min] 
Ingress duration, Ti ng [min] 
R p [i? ffi ] 
M p [M©] 



9.2869944±0.0000088 a 
4959.70597±0.00036 
0.02379±0.00012 



21.59; 
89.94" 



-0.15 
-0.37 
-0.06 
-0.44 



0.466±0.061 
0.155±0.097 

u - uz -0.02 

201.91±0.47 

4 7n+ - 18 

2.209±0.048 
<20.3 b 



Kepler- 19c 


Period [days] 


<160 c 


M p [M Jupitcr ] 


<6.0 C 



a The period and T c values for Kepler-19b are deter- 
mined from a linear fit to the transit times. 

b The upper limit on the mass of Kepler-19b is de- 
termined from the radial velocity analysis in Section 
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4.2. 

c The upper limit on the mass and period of Kepler- 
19c is described in the text of Section 6.1.2. 
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4. Follow-up Observations 

4.1. Reconnaissance Spectroscopy 

We gathered reconnaissance spectra of Kepler-19 on UT 2009 August 05 (orbital phase 0.596), 
2009 August 29 (phase 0.161) and 2010 October 1 (phase 0.017) with the Tull coude spectrograph 
of the McDonald Observatory 2.7m telescope. We used these high resolution (R=60,000) spectra 
to verify the Kepler Input Catalog stellar classification and to search for evidence of any secondary 
stellar spectrum or binary orbital motio n. We cross-co r relate d the spectra against a library of 



synthetic stellar spectra as described by iBatalha et al.l (|201ll ). The spectra did not show any 



evidence of a secondary spectrum, and the absolute radial velocities, which cover both orbital 
quadratures, agree at the 0.75 km/s level. These spectra gave the best match to the synthetic 
templates with T e g- = 5750 K, log(g)= 4.5, and v sin i = 2 km s _1 , for an assumed solar metallicity. 
The height of the cross-correlation peaks was 0.98 for all of the spectra, indicating an excellent fit 
to the stellar template spectra. 



4.2. High-resolution Spectroscopy 



Between 2009 October 29 and 2011 June 10 we gath ered 8 high-resolution spectra of Kepler- 
19 with the Keck HIRES spectrometer (jVogt et al.lll994l ). With these spectra, we conducted an 
analysis to determine the stellar parameters. We compared a high-resolution template spectrum 
to stellar models, generated with the spectral syn thesis package Spectroscopy Made Easy (SME; 
Valenti &: Piskunovl Il996l : IValenti &: Fischerl 120051 ) . We determine the effective temperature T e fj 
of the host star of 5541 ± 60 K, a metallicity [Fe/H] of -0.13 ± 0.06 dex, a log(g) of 4.59 ± 
0.1 dex, and a v sin i < 2 km s . We comment briefly here on the stellar activity. We find a 
value of the ratio of emission from the Ca II H and K lines to the total bolometric emission of 
log(ilffy)=-4.95±0.05. T he R H k value is derived from a Mt. Wilson style S-value of 0.174±0.007 
([Isaacson Fischerll2010i ). The log(RjjK) value is low for main sequence stars of this temperature 
and is consistent with the slow stellar rotation we infer from the measured v sin i of <2 km s -1 . If 
we assume rigid body rotation of the star and a stellar spin axis aligned with the orbital spin axis 
of the planet, we find a lower limit on t he stellar rotation period of 22 days. The rotation period 
derived from the Rhk value is 32 days (jNoves et al.lll984l ) which along with the lack of emission 
in the core of the Ca II H and K lines leads us to conclude that the star is relatively inactive. 

We further used the spectra to derive estimates of the stell ar radial velocity. The spectra 
were gathered with the same configuration of HIRES, described in lMarcv et al.l (|2008l ). which was 
demonstrated to yield typical precisions of 1.0 to 1.5 m s _1 on nearby FGK stars. This method 
relies on the use of an iodine cell placed in front of the beam, which superimposes the iodine 
spectrum on the stellar spectrum with an identical instrumental profile. For each 100 pixel section 
of the spectrum, the iodine and stellar spectral lines are fit simultaneously. For the set of these 
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observations, this treatment yields a typical internal error estimate on individual radial velocities of 
1.5 m s . We note also the use of the "C2 dec ker" entrance apertur e for all of these observations, 



which technique is described in greater detail in lBatalha et aJj (|201ll ). and enables sky subtraction 



(as compared to the "B5 decker", for which sky subtraction is not possible). We list the measured 
radial velocities, with associated error bars (excluding stellar jitter), in Table [3j We gathered 
twelve additional observations prior to the ones listed in Table [3] but these were observed with the 
B5 configuration, and had a much higher scatter (15 m s _1 as opposed to 4 m s -1 ). For this reason, 
we excluded them from the analysis. 

We determine the upper mass limit on Kepler- 19b from the ra dial velocities as follows. We 



employ the Bayesian MCMC technique described in iGregorvl (|2007l ) to fit a radial velocity model 
to the observations. The free parameters in the model are the semi-amplitude velocity K of the 
star, the zero-point velocity 7, the eccentricity e of the planetary orbit, the argument of perihelion 
u, and a stellar jitter term. The orbital period and transit epoch are also free parameters, however 
the precision of the Gaussian priors we place on them from the light curve analysis (see Table [1]) 
effectively fixes their values. Additionally, we fix the inclination of the orbit to the value measured 
from the light curve of 89.94°. We first fit a model, assuming a circular orbit. In this case, we infer 
a stellar jitter contribution of 4.1 ± 1.7 m s , and 1, 2 and 3a upper limits on the semi-amplitude 
of 1.4, 3.3, and 4.9 m s _1 , respectively (these values are derived from integrating over the posterior 
distribution of semi-amplitude until 68%, 95%, and 99.73% of the area is enclosed). Applying the 
semi-amplitude toward a mass upper limit, we determine a 3a upper limit on the mass of < 15.2 
M®, if the orbit is circular (these upper limits are 4.3 and 10.2 M® at 1 and 2a confidence). The 
most likely fit (depicted as the solid line in Figure [6|) has an amplitude of 0.5 m s _1 (corresponding 
to a mass of 1.6 M®), but this value is well below our detectability threshold. For comparison, 
we also show representative circular orbit at the 3a upper limit for semi-amplitude (dashed line in 
Figure [6|) . 

We also address the possibility of a non-zero eccentricity. If the orbit were significantly ec- 
centric, then the transit duration of Kepler-19b would deviate from the predicted duration for the 
edge-on circular orbit scenario, unless the argument of perihelion conspired to mimic the circular 
transit duration. As eccentricity increases, there is an increasingly narrow range of oj that matches 
the transit duration of a circular orbit. As we demonstrate in Section 3, the assumption of a circu- 
lar orbit is consistent with the well-constrained low impact parameter (high inclination) measured 
from the the Kepler light curve, as illustrated by the short ingress and egress times. In addition 
the log(g) inferred from the light curve analysis assuming zero eccentricity agrees with the inde- 
pendent spectroscopically determined value. We note that the planet is too small to constrain the 
eccentricity from secondary eclipse observations; we discuss this possibility further in Section 6.3. 
Therefore, we elected to impose a prior on e and u from our knowledge of the transit duration, 
as follows. The ratio r of the transit duration for an eccentric orbit to the transit duration for a 
circular or bit can be ap proximated by the following expression of the eccentricity and argument of 



perihelion (lBurkell2008h : 
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r = (1 - e 2 ) 1 /2/( 1 + e cos ( w - vr/2)). (3) 

For each element of the MCMC chain, for which we now vary e and uj, we evaluate the transit 
duration ratio r. We then assign a flat prior on r, which is equal to one for 0.7< r <1.3 (that is, 
the transit duration is within 30% of circular) and zero otherwise. Applying this transit duration 
prior during the radial velocity parameter estimates, we find a negligibly smaller estimate for the 
stellar jitter than the circular orbit case and semi-amplitude upper limits (again, with 1, 2, and 
3a confidence respectively) of 2.0, 9.2, and 23.5 m s _1 . These semi-amplitudes are associated with 
mass upper limits of 5.5, 20.3, and 50.3 M®, with the same stated confidences. Given the radius 
value for the planet determined from the Kepler photometry and the 2a upper mass limit of 20.3 
M®, we find an upper limit on the density of Kepler-19b of 10.4 g cm -3 . We comment further on 
the possible composition of the planet, given these upper limits, in Section 6.4. 

10 F ' ' ' ' ' ' 1 
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Fig. 6. — Measured radial velocities of Kepler-19, as a function of phase (assuming the orbital 
period and epoch of Kepler-19b, as stated in Table [1]). We have depicted two radial velocity models 
with zero eccentricity, the first (dashed) corresponding to a planetary mass at the 3a upper limit 
of 15.2 M®, and the second (solid) corresponding to the most likely amplitude of 0.5 m s _1 (or 1.6 
M®). The error bar depicted include the effects of stellar jitter, which we conclude are near 4 m 
s- 1 . 



4.3. Adaptive Optics Imaging 



We gathered adaptive optics i mages in J band of Ke pler-19 on UT 24 September 2009, using 



the PHARO near-i nfrared camera (Havward et al.l l200ll ) on the the Hale 200 inch telescope on 
Mt. Palomar, CA. I Troy et al.l ( 20001 ) give a complete description of the Palomar adaptive optics 
system. We employe d a dither pattern for these observations similar to the technique described in 



Batalha et al.l (|201ll ). although we used an exposure time of 4.25 seconds. In Figure [71 we show 



the local neighborhood of Kepler-19 within 10 arcseconds. 
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Fig. 7. — J-band adaptive optics image of the neighborhood of Kepler-19, within 10 arcseconds 
(left) and 1 arcsecond (right). 



We assess our s ensitivity to additional sources using a similar procedure to that described by 
Batalha et al.l (|201ll ). We inject fake sources near the target star at random position angles, using 
steps in magnitude of 0.5 mag and varying the distance from the target star in increments of 1.0 
FWHM of the point-sp read function (PSF). We then attempt to identify the injected sources with 
the DAOPhot routine (|Stetsonlll987l ) and also by eye, and set our sensitivity limit, as a function of 
distance, at the magnitude where we are able to recover the injected sources. The limit in Am as a 
function of distance from the target star is shown in Figure [HJ We then convert the Am sensitivity 
limit in J band to a limit in Kepler magnitudes, by assuming a nominal Kepler magnitude- J color 
(using the value derived from a magnitude-limited sample of Kp- J=1.28±0.52 mag). We do not 
detect any additional sources within our sensitivity limits in the neighborhood of Kepler-19. 



4.4. Speckle Imaging 



We gathered speckle images of Kepler-19 using filters in both R and V band on UT 18 June 
2010, using the D ifferential Speckle Survey Instrument located at the WIYN teles cope (DSSI, 



Hor ch et al.ll2009i ) A detailed discussion of the recent upgrades to DSSI is presented in lHorch et al 



(|2010l ). and summary of the speckle ima ging survey of Kepler candidates, and those reduction 
procedures, is given by Howell et al.l (j201ll ). We assess our sensitivity to the presence of additional 
stars near the Kepler target star as a function of angular distance. For concentric rings of varying 
radius, centered on the target star, we determine the magnitude difference between the target star 
itself and the local extrema of the sky background. Figure [9] shows the results of this procedure 
in both R and V band. We find that we would have detected a companion at a distance of 0.2" 
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Radial Distance [arcsec] 



Fig. 8. — The sensitivity limits to additional point sources in the neighborhood of Kepler-19b as a 
function of radial distance from the primary target. The filled circles represent the J-band limits 
and each point represents a step in FWHM away from the primary target centroid peak. The 
dashed line underneath represents the J-band limits converted to Kepler magnitude limits if a star 
were to have a nominal Kp-J color, as described in the text. 



with a difference in magnitude smaller than Am=3.48 in R band, and a companion at a distance 
of 0.25" with Am <1.65 in V band. Here again, as in the adaptive optics images, we detect no 
additional sources nearby to the Kepler target star. 



5. Planetary Validation of Kepler-19b 

5.1. Photocenter Tests 

We use two methods to search for false positives due to background eclipsing binaries, based on 
examination of the pixels in the aperture of Kepler-19: direct measurement of the source location 
via difference images, and inference of the source location from photocenter motion associated with 
the transits. We employ two methods because of their different vulnerabilities to systematic bias; 
when the methods agree, we have increased confidence in their result. 



Difference image analysis ([Torres et al. 



2011 



takes the difference between average in-transit 
pixel images and a verage out-of-transit images. A fit of the Kepler pixel response function (PRF; 
Brvson et al.ll2010i ) to both the difference and out-of-transit images directly provides the location 
of the transit signal relative to the host star. We measure difference images separately in each 
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Fig. 9. — Top panel: R band speckle sensitivity curve of Kepler-19. The magnitude difference be- 
tween the target star and local extrema in the background are denoted by squares (local maxima) 
and points (local minima). The solid line denotes a flux that is 5a brighter than the mean back- 
ground brightness (where a is the standard deviation of the extrema in the background), where we 
could confidently detected an additional source. No companions are detected within 1.8 arcseconds 
of the target star to a depth of 4 magnitudes. Bottom panel: V band speckle sensitivity curve of 
Kepler-19. No companions are detected within 1.8 arcseconds of the target star to a depth of 2 
magnitudes. 
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quarter, and estimate the transit source location as the robust uncertainty- weighted average of the 
quarterly results. 

We measure photocenter motion by computing the flux-weighted centroid of the pixels in the 
optimal aperture, plus a one-pixel halo in every cadence, generating a centroid time series for row 
and column. We fit the modeled transit to the whitened centroid time series transformed into sky 
coordinates. We perform a single fit for all quarters, and then infer the sou rce location by scalin g 



the difference of these two centroids by the inverse of the flux as described in I Jenkins et al.l (|2010al ) 



The source as determined by the difference image method is offset from the nominal location 
of Kepler-19, as given in the Kepler Input Catalog, by 0.09 ± 0.11 arcsec = 0.80a. The source as 
determined by the flux-weighted centroid method is offset from Kepler-19 by 0.10 ± 0.11 arcsec 
= 0.88o\ The location of the offsets is shown for both methods in Figure [TUJ Both methods 
show that the observed centroid location is consistent with the transit occurring at the location of 
Kepler-19. 




Fig. 10. — Quarterly and average reconstructed transit source locations relative to Kepler-19. Left: 
The green crosses show the individual quarter measurements using the difference image technique, 
and the magenta cross shows the uncertainty-weighted robust average of the quarterly results. 
Right: The magenta cross shows the transit source location reconstructed from the multi-quarter 
fit of the transit signal to the centroid motion. The length of the crosses show the la uncertainty 
of each measurement in RA and Dec. The circles show the 3a circle around the average source 
location. The location of Kepler-19 is shown by the red asterisk along with its Kepler ID and 
Kepler magnitude. 
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5.2. Spitzer Observations 



Warm Spitzer observations in the near-infrared can al so prove useful toward validating Kepler 

Unless a putative blend scenario is 



candidates, as shown for Kepler-lOc (jFressin et al 



2011 



comprised of stars of nearly identical color, the transit depth in a blend scenario will depend upon 
the wavelength at which it is observed. Conversely, an authentic transiting planet will produce an 
achromatic transit depth. 



We gathered observations using the Infrared Array Camera (IRAC) (jFazio et al.l 12004 ) on 
Warm Spitzer at 4.5 fim of two consecutive transits of Kepler-19: one on UT 29 June 2010, and 
one on UT 9 July 2010. Both observations span 8 hours, centered on the 3.5-hour-long transit. 
We gathered the observations using the full- array mo de of IRAC , with an integration time of 26.8 
s/image. We employed the techniques described in lAgol et al.l (|2010l ) for the treatment of the 
images before photometry. We first converted the Basic Calibrated Data products from the Spitzer 
IRAC pipeline (which applies corrections for dark current, flat field variations, and detector non- 
linearity) from mega-Janskys per steradian to data number per second, using 0.1469 MJy-sr -1 per 
DN s _1 , and then to electrons per second, using the gain of 3.71 e DINT 1 . We identified cosmic 
rays by performing a pixel-by-pixel median filter, using a window of 10 frames. We replace pixels 
that are > 4<r outliers within this window with the running median value. We also corrected for 
a striping artifact in some of the Warm Spitzer images, which occurred consistently in the same 
set of columns, by taking the median of the pixel values in the affected columns (using only rows 
without an overlying star) and normalizing this value to the median value of neighboring columns. 

We estimate the position of the star on the array using two techniques. First, we employed 
a flux-weighted sum of the position within a circular aperture of 3 pixels (we tested whether the 
size of this aperture made a difference by increasing the size to 4 pixels and repeating the analysis: 
we found that the position estimates were nearly identical). Additionally, we fit a Gaussian to the 
core of the PSF using the IDL routine GCNTRD (again using apertures of 3 and 4 pixels, and finding 
no significant difference between them). We then performed aperture photometry on the images, 
using both estimates for the position and variable aperture sizes between 2.1 and 4.0 pixels, in 
increments of 0.1 pixels up to 2.7 pixels, and then at 3.0 and 4.0 pixels. We decided to use the 
position estimates using a flux-weighted sum at an aperture of 2.6 pixels, which minimized the 
out-of-transit RMS. Alternatively, using the positions derived with GCNTRD, using a slightly smaller 
aperture, or using a slightly larger aperture, changed the RMS by only a few percent at the most. 

We rem ove the effect of the IRAC intrapixel sensitivi ty variations, or the "pixel-phase" ef- 
fect (see eg. ICharbonneau et al.ll2005l : iKnutson et al.l 120081 ) using two techniques. With the first 
technique, we assume a polynomial functional form for the intrapixel sensitivity (which depends 
upon the x and y position of the star on the array). We denote the transit light curve / (which 
depends upon time), a nd we hold all paramet ers constant except for the transit depth. We use the 
light curve software of Mandel fc Agoll ( 2002 ) to generate the transit models. The model for the 
measured brightness f'(x,y) is given by: 
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/' = f(t, R p /R+) • (h + b 2 (x -x) + b 3 (x- x) 2 + b 4 (y -y) + b 5 (y- y) 2 ), (4) 

where we include all of the observations (both in- and out-of-transit) to fit the polynomial 
coefficients and the transit depth simultaneously. 

We included cross-terms in x and y, as well as higher order terms, but found that did not 
substantially decrease the RMS error of the out-of-transit residuals after the sensitivity function 
is divided from the flux. We fit for the polynomial coefficients b\ through 65 using a Levenberg- 
Marquardt x 2 minimization. We also performed an MCMC analysis to fit the polynomial co- 
efficients to determine whether fitting for the transit depth was degenerate with any other free 
parameters, and determined that about 20% of the error in the best-fit transit depth is due to a 
degeneracy with the strongest polynomial dependence of the intrapixel sensitivity, which is the lin- 
ear coefficient in y. However, the Spitzer light curve contains significant correlated noise even after 
the best intrapixel sensitivity model is removed. We incorporate the e ffect of rema i ning correlated 



noise with a residual permutation analysis of the errors as described bv lWinn et al.1 (|2008l ). wherein 
we find the best-fit model /' to the light curve as given by Equation [H subtract this model from 
the light curve, shift the residuals by one step in time, add the same model back to the residuals, 
and refit the depth and pixel sensitivity coefficients. We wrap residuals from the end of the light 
curve to the beginning, and in this way we cycle through every residual permutation of the data. 
We determine the best value from the median of this distribution, and estimate the error from the 
closest 68% of values to the median. We gathered 4.5 hours of observations outside transit, which 
is sufficient to sample the systematics on the same timescale as the 3.5 hour transit. Using the 
residual permutation method on the light curve treated with a polynomial, we find a best-fit transit 
Rp/R* to be 0.0226 ± 0.0039 for the visit on 29 June, which is consistent with the best solution 
using MCMC, although the error bars are inflated by 40% when compared to the MCMC error 
bars. For the visit on 9 July, we find 0.0280 ± 0.0027 with the rosary-beading analysis; these error 
bars are 20% larger than the corresponding MCMC error bars. The larger error on the transit 
depth measured on the first visit is attributable to the larger area on the pixel over which the star 
wanders during the observations: the smaller this area, the better we are able to fit the polynomial 
sensitivity model. While the extent of the pointing oscillations in the x direction are comparable 
between the two visits (0.1 pixels), they differ substantially in the y direction. The star wanders 
0.15 pixels in y on the 29 June visit and 0.08 pixels in y on the 9 July visit. 



We also treated the light curve with the weighted sensitivity function used in iBallard et al 



(|2010bl ). which proved in that work to produce a time series with lower RMS residuals. For this 
procedure, we do not assume any a priori functional form for the intrapixel sensitivity; rather, 
we perform a weighted sum over neighboring points for each flux measurement, and use this sum 
to correct each flux measurement individually. In this way, we build up a map of the intrapixel 
sensitivity point-by-point. We use same widths, cr x =0.01 7 pixels and cx^O.OO^ S, for the weighting 



function (which is a Gaussian in x and y) as we used in lBallard et al.l (|2010al ). We therefore have 
only one free parameter in this case, which is R p /R+. We correct each observation using all other 
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flux measurements, but we did not bin the data, as was done in lBallard et al.l (|2010al ). With this 
treatment, we also use the residual permutation method to fit the transit depth at each residual 
permutation. We find i? p /^=0.0242±0.0032 for the 29 June visit and i? p /i2*=0.0233±0.0033 for 
the 9 July visit. These errors are larger by 20% and 30%, respectively, as compared to the MCMC- 
derived errors. The out-of-transit RMS of the light curve is slightly lower in both cases using the 
weighted sensitivity function treatment, so in this case we defer to the weighting- function-derived 
values for R p /R+. Combining the two measurements, we find a radius ratio R p /R± at 4.5 /mi of 
0.0238±0.0023, which translated to a transit depth of 5471}^ ppm. This value is in excellent 
agreement with the depth in the Kepler bandpass of 567±6 ppm (corrected for limb-darkening). 
In Figure fTTj we show the combined and binned Spitzer light curve, with the best-fit transit model 
derived from the Spitzer observations and the best-fit Kepler transit model (corrected for limb- 
darkening) overplotted. We comment further on the types of blends we rule out with Spitzer (and 
their complementarity with blends ruled out by BLENDER) in the following section. 



5.3. BLENDER Analysis 



In the absence of a radial velocity confirmation and mass measurement of the planet Kepler- 
19b, we instead investigate the likelihood that the transit signal is a false positive. Possible false 
positive scenarios involve another eclipsing system lying within the same photometric aperture as 
the Kepler target star. This binary system could comprise two stars, or a star and a gas giant 
planet, and could be physically associated or unassociated with the host star (in the foreground or 
background). In a false positive scenario, the presence of the Kepler target star dilutes the depth of 
the transit to appear planetary (or attributable to a smaller planet, if the binar y system comprises 
a sta r and a Jupiter-size planet). We employ the BLENDER software package ([Torres et al.l 120041 . 



2011 



), which produces synthetic light curves corresponding to eclipsing binary blend scenarios and 
attempts to replicate the detailed shape of the Kepler transit light curve. The BLENDER techniqu e 



has been ap 
Kepler- 1 If ( 



plied previously towa rd validating three Kepler planets: Ke pler-9d ([Torres et al 



Lissauer et al 



2011 



2011 



2011 



Model blend light curves 



and Kepler- 10c (jFressin et al. 
are compared with the Kepler photometry in a x 2 sense, with models considered poor fits accord- 
ingly deemed unlikely to explain the transit. By exploring the parameter space of mass, impact 
parameter, orbital eccentricity, and distance from the host star, BLENDER amasses knowledge of 
the range of possible blends that are consistent with the shape of the transit, and the range of 
blends that are inconsistent. The BLENDER nomenclature defines the objects within the binary to 
be the "secondary" and "tertiary" , while the Kepler target star is defined to be the "primary" . In 
the case of hierarchical triples, BLENDER uses the best isochrone parameters for the Kepler target 
star (derived from SME) as input constraints to the secondary and tertiary stars: these stars are 
assumed to have formed concurrently and are therefore assumed to be the same age. When the 
secondary and tertiary are physically unassociated wit h the Kepler target star, BLENDER assumes 
an age of 3 Gyr (a representative age for the field, per lTorres et al.ll201ll ) and a solar metallicity 
to model the putative binary system. 
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Fig. 11. — Both transits of KOI-084 gathered with Warm Spitzer at 4.5 ^m. The top panels show 
the raw flux, binned by a factor of 4, with the intrapixel sensitivity variation (obtained with the 
weighted sensitivity function, as described in the text) overplotted in red. The middle panels show 
the individual transits with this intrapixel sensitivity removed, binned by a factor of 16, with the 
best models overplotted. The bottom panel shows the combined transit light curves gathered with 
Spitzer. The best-fit transit model with depth derived from the Spitzer observations is shown in 
red, while the Kepler transit model (corrected for limb darkening) is shown in green. The Spitzer 
and Kepler transit depths are in excellent agreement. 



In order to validate Kepler-19b as a planet, we evaluate the probability of false positive sce- 
narios allowable by BLENDER and compare these probabilities to that of the authentic 2.2 R§ planet 
hypothesis. First, we address the probability of a physically unassociated binary in the foreground 
or background of the target star. In the case of Kepler-19b, all false positive scenarios consisting of 
a background or foreground eclipsing binary (comprising two stars) are ruled out at the 10a level 
from the shape of the transit alone. That is, the best planet model furnishes a solution that is a 
10fT improvement over the best blend model in this case. This result is attributable to the sharp 
ingress and egress of the transit light curve, which is not well reproduced by blend models involving 
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a binary system consisting of two stars. We therefore confidently conclude that this scenario cannot 
replicate the Kepler transit signal. 

For scenarios consisting of a foreground/background star orbited by a larger (Jupiter-size) 
planet, there exists a region of parameter space in which an eclipsing binary model provides a 
comparably good fit, as compared to the single star and planet model. The BLENDER constraints 
are represented by contours of equal goodness of fit in Figure [12j The 3c contour is shown with 
a heavy white line, and blend scenarios under this curve are considered viable. In additional to 
the goodness-of-fit of the blend models to the Kepler light curve, there are regions of parameter 
space that are disallowed by the color of the star (as measured by 2MASS and in Sloan r band), 
as well as the spectrum (in which a secondary star within a certain magnitude range would be 
apparent). These constraints are depicted in Figure [T2l as blue cross-hatching (within which region 
blends are disallowed) and a solid green line (below which an additional star would have appeared 
in the spectrum). We comment briefly on this spectroscopy constraint, which we measured by 
injecting additional stellar spectra (in this case, solar-type) into the spectrum of Kepler-19 at 
varying brightnesses and relative velocities, and then determining the limits on detectability via a 
cross-correlation of the spectrum with a template. We determine that any star within 10% the flux 
of Kepler-19 would be detectable in the cross-correlation function down to relative velocities of 5 
km s _1 , which translates conservatively to a Am < 2 constraint. A velocity variation of <5 km s _1 
would be unlikely for a random unassociated background star, and we comment on the hierarchical 
triple case further below. Furthermore, the angular separation of this star and planet system must 
also be sufficiently small as to be undetectable by adaptive optics imaging for which limits are 
given in Section 14.31 While an unassociated binary might be resolvable by adaptive optics, for a 
hierarchical triple the possibility of an unresolved companion remains. 

We proceed to evaluate the frequency of the remaining allowable blends as follows: We estimate 
a priori the frequency of stars in the background or foreground of the target star. We evaluate a 
theoreti cal number d e nsity of neighboring stars, per square degree, using Galactic structure models 
given bv lRobin et al.1 (|2003l ). We record this number density in half-magnitude bins (shown in Table 
2]), and show both the number density of stars and their allowed mass values in Columns 3 and 4, 
based on the constraints constraints imposed by the BLENDER contours, additional color constraints, 
and brightness constraints. For the magnitude bins in which no blend furnishes a solution within 
3a of the planetary model, we have left these columns blank. The maximum angular separation 
at which these stars might have remained undetected in adaptive optics imaging (the limits are 
provided in Section [4.3p is listed in Column 5. The number of stars in each magnitude bin is then 
listed in Column 6. In order to evaluate the tr ansiting planet prior, we rely on the reported Kepler 
planet candidate sample to date, presented in iBorucki et al.l (|201ll ). While the majority of these 
candidates have not ye t been confirmed to be pl anets, the false positive rate is expected to be quite 
small (as reported bv iMorton fc Johnson! 120 111 ) and so will no t subs tantially change our results. 



We also assume that the sample presented in the IBorucki et al.l (|201ll ) candidate list is complete. 



The second feasible blend scenario is an additional star and transiting Jupiter-size planet, 
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Secondary mass (M ) 



Fig. 12. — BLENDER \ 2 goodness-of-fit contours corresponding to blend models with back- 
ground/foreground (physically unassociated) secondary star and planetary tertiary, as a function 
of distance modulus and mass of the secondary. Viable blend scenarios lie below the 3a contour, 
depicted with a heavy white line. The blue cross-hatching shows the region of parameter space 
for which the blend model is the wrong color to be consistent with the measured 2MASS and 
Sloan r colors of the star, while models that lie below the solid green line have a small enough 
magnitude difference that the secondary would have been detected in the spectrum of the star 
(Am < 2.0). The dashed green line shows location of the Am contrast limit corresponding to the 
faintest allowable blend. The remaining allowed parameter space is denoted "Allowed Region." 

which are physically associated with the Kepler target star. We present the results of BLENDER for 
this case in Figure [131 While there exists a range of hierarchical triples whose light curve shape is 
consistent with the Kepler transit (depicted by the x 2 contours) , these are all ruled out by either 
the color constraint on blends (shown in blue cross-hatching) or the brightness constraint on blends 
from the spectrum (shown in green cross-hatching). There exists the remaining possibility of a 
true twin to the target star: a star that has an identical color, and whose position is either at a 
distance > 20 AU (at which position the predicted radial velocity is equal to 5 km s _1 ) or whose 
tangential velocity is <5 km s _1 relative to the target star during the time of our observations. 
Even if this scenario were to be true, the planet's inferred radius would only be greater by a factor 
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of \/2. Therefore, we conclude that the only possible blends belong to the unassociated planet and 
star scenario. 




0.2 0.4 0.6 0.8 

Secondary moss (M ) 

Fig. 13. — BLENDER x 2 goodness-of-fit contours corresponding to hierarchical triple blend models 
(with physically associated secondary star and planetary tertiary). As in Figure [T21 viable blend 
scenarios lie within the 3a contour, depicted with a heavy white line. The color-coding is similar 
to the previous figure, with blends of the wrong color depicted in blue cross-hatching, and blends 
that are ruled out spectroscopically depicted in green cross-hatching. 



Combining the probabilities associated with all background or foreground star and planet pairs, 
we find a total blend frequency of 1.08xl0~ 7 . This frequency corresponds to the likelihood of a 
blend that is capable of producing a transit light curve that is no worse (< 3a) than the best model 
corresponding to a transiting planet around a single star within the aperture. 

Next, we estimate the a priori frequency of a true planet with the characteristics implied by 
the Kepler transit light curve. Using the planetary radius range of 2.209±0.048 i?®, we identify 



119 planet candidates from the lBorucki et al.l (|201ll ) catalog with sizes that are within 3a of this 
measured value. The planet prior is equal to 7.6x 10~ 4 (119 divided by the total number of Kepler 
targets, 156,453), which is more than 3 orders of magnitudes larger than the blend frequency. We 
therefore find that the planetary scenario corresponding to a 2.2 Rq planet is 7000 times as likely 
as the blend scenario, and conclude with very high confidence that the transit signal is due to a 
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planet, Kepler-19b. 

For comparison, the constraint from Spitzer in this case provides an independent means of 
ruling out a subset of blends which are also ruled out by BLENDER. If we impose the constraint 
that a putative additional star may not produce a transit depth at 4.5 //m which is 3a deeper than 
measured, such a star cannot be less massive than 0.7 M & (otherwise, the additional star would be 
so red as to produce a significantly deeper transit depth in the near- infrared) . As shown in Figures 
[12] and [131 BLENDER independently rules out blends consisting of a star in this mass range from the 
shape of the transit light curve. 
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Table 3. Relative Radial Velocities for Kepler-19 



HJD 


RV 


Unc. 


-2450000 


m s _1 


m s _1 


5134.805 


0.9 


2.1 


5320.113 


0.4 


1.7 


5402.956 


5.9 


1.5 


5407.827 


-6.5 


1.4 


5412.981 


-0.2 


1.6 


5435.879 


5.5 


1.3 


5723.074 


2.0 


1.5 


5723.950 


4.0 


1.5 



Note. - Uncertainties do 
not include stellar jitter, which 
is likely to be near 4 m s _1 . 
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Table 4. Blend frequency estimate for KOI-084.01. 



Blends Involving Planetary Tertiaries 



Kp Range 


AKp 


Stellar 


Stellar Density 


Pmax 


Stars 


Transiting Planets 


(mag) 


(mag) 


Mass Range 


(per sq. deg) 


(") 


(xl0~ 6 ) 


0.36-2.00 R 3up , /pi an et = 0.20% 






(M e ) 








(xlO- 6 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


11.9-12.4 


0.5 












12.4-12.9 


1.0 












12.9-13.4 


1.5 












13.4-13.9 


2.0 












13.9-14.4 


2.5 


0.88-1.40 


444 


0.22 


5.21 


0.011 


14.4-14.9 


3.0 


0.91-1.40 


505 


0.29 


10.3 


0.021 


14.9-15.4 


3.5 


0.95-1.40 


436 


0.38 


15.3 


0.031 


15.4-15.9 


4.0 


1.00-1.30 


327 


0.53 


22.3 


0.045 


15.9-16.4 


4.5 












16.4-16.9 


5.0 












16.9-17.4 


5.5 












17.4-17.9 


6.0 












17.9-18.4 


6.5 












18.4-18.9 


7.0 












18.9-19.4 


7.5 












19.4-19.9 


8.0 












19.9-20.4 


8.5 












Totals 






1712 




53.11 


0.108 



Total frequency (BF) = (0.108) X 10~ 6 = 1.08 X 10~ 7 



Note. - 



Magnitude bins with no entries correspond to brightness ranges in which BLENDER excludes all blends. 
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6. Discussion and Conclusions 

6.1. Interpretation of Transit Timing Variations 

Our analysis of the transit timing variation of Kepler-19b comprises two sections. In Sec- 
tion ETTJ w e argue that there must exist a second planet in the Kepler-19 system, since the TTVs 
cannot originate from astrophysical effects or a stellar mass perturber. Then, in Section [6,1.21 we 
describe the dynamical properties of planetary perturbers that could account for the observed TTV 
pattern of planet b. 



6.1.1. Demonstration of Planethood of Perturber 



In Figure 4 we presented a strongly detected and nearly sinusoidal variation (period Pttv = 
316 days and amplitude A t t v = 5 min) in the times of transits of planet b. Here we discuss four 
potential interpretations of the signal which do not invoke a perturbing planet, and demonstrate 
that the signal cannot originate from these scenarios. The first two scenarios consider a system 
of only the star and the transiting planet, while the latter two allow for the presence of a third, 
non-planetary, body. Because these scenarios are disallowed, the only viable alternative is that the 
signal is a dynamical variation caused by a second planet (planet c), which we discuss in the next 
subsection. 

First, we consider the possibility that the signal is due to stellar activity. The most plau- 
sible c andidate for dynamical interaction with the star causing the TTV signal is the lApplegate 
(Il992l ) effect from the eclipsing binary literature, which IWatson &; Marshl (|2010l ) recently applied 
to exoplanets. In this mechanism, the star undergoes a magnetic cycle, which varies the rotational 
bulge's gravitational pull on the planet, slightly varying its orbital period. For this effect to produce 
the observed TTV signal, the magnetic cycle would need to have an exceptionally short period of 
316 days. This is problematic because the stellar activity in the spectra is low, as we describe in 
Section 4.2, suggesting a magnetically in active star with a l o ng sp in period. For magnetic dynamos 
typical of solar-type stars like Kepler-19, IWatson &: Mars h (2010) calculate TTV variations of less 
than 1 second over timescales of several years, much too small to explain our data. Apart from this 
dynamical interaction, activity could cause apparent TTV via the p l anet t ransiting over starspots 
(|Silva- Valid 120081 ; lAlonso et al.ll2008l ). For instance, iKnutson et al.l (|201ll ) recently found transit 
timing deviations for GJ 436b that are greater in optical photometry, where spots are pronounced, 
than in near-infrared photometry, where spots are relatively muted. However, in our case, even 
with excellent Kepler photometry, no spots are detected, either in out-of-transit stellar modulation 
or in excess residuals during transit due to spot crossings. 

Second, the signal c ould be due to rotation of the planetary orbit's apsidal line about the star 
(|Hevl fe Gladmanl 120071 ). The eccentricity need not be large (e b ~ 2irA t t v /P b = 0.0023), but the 
apse must be precessing very quickly to be consistent with the TTV period (Pttv = 316 days). 
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Using expressions for realistic apsidal motion rates (|Fabrvcky||201Cll . section 3.1.1), we demonstrate 
that the periods associated with possible precession mechanisms are too long by several orders of 
magnitude, as follows. General relativistic precession has a period of 7.9 x 10 yr. A star made 
oblate by rotation, alternatively, would cause precession wit h a period of 7x 10 6 yr x( P rot /10 days) 2 , 
if the star has an apsidal motion constant of &l/2 = 0.02 (jClaret &: Gimenedll992l ). Finally tidal 
distortion of the plan et would cause precess ion with a period of ~ 10 8 yr, if the planet has a Love 
number of = 0.3 (jMardling fc Linl 12004 ). We conclude that general relativity dominates the 
putative prec ession rate, but it is inc onsistent by orders of magnitude with P ttv and is probably 
undetectable (jRagozzine fc W olf 2009). 



Third, the signal could originate from light time delay, owing to reflex motion of the system 
as it is orbited by a third body with an orbital period of P 2 = Pttv — 316 days. It is conceivable 
that this putative body moves the bary center of the Kepler- 19 / Kepler- 19 b system by ±5 light 
minutes, causing a time-variable light-time delay (e.g.. llrwinlll952l : ISvbilski et al.ll2010l ). In creating 
such a large displacement in the radial direction, a radial velocity signal - its derivative - would 
also be created. The semi-amplitude would be 27r(0.6AU)/(316days) ~ 17 km s _1 , ruled out by 
orders of magnitude given our radial velocity measurements. Moreover, to induce this motion, the 
additional body must have a mass of at least I.OMq, and would likely impart a second set of lines 
on the spectra, which are not detected. 

Finally, the signal could be dynamical, owing to perturbations from another body in the system 
such as a second star or a brown dwarf. The orbital period of such a body cannot be too long; 
otherwise, the gravitational potential it induces on planet b would result i n a period longer than the 
observed Pttv In fact, the longest period the body could have is 2P ttv (jBorkovits et al.ll2003l ). as 
the dominant part of the perturbation of a distant body is its "tidal" term, which has a frequency of 
twice the body's orbital frequencjQ. The lack of large radial velocity variation means that additional 
bodies must have masses in the planetary regime (a case handled in the next subsection) on nearly 
any orbit with P2 < IPttv An exception is for nearly face-on orbits with a small component of 
radial motion. To our knowledge, there is only one such configuration that could explain the timing 
data, as follows. A circula r orbit of period 2Pu v , mutually perpendicular to the orbit of planet b, 
causes a timing signature (jBorkovits et al.ll2003l . eq. 46) : 



o-c b 



M 2 



Pi 



sin 2/2 



(5) 



3tt M± + M b + M 2 P 2 

where M 2 is the perturber's mass and f 2 is its true anomaly measured with respect to the plane 
of planet b's orbit. The perturbing orbit must be nearly circular in the case of Kepler-19, as we 
have measured two full O — C cycles and found them to be nearly identical, rather than different 
ei ther in amplitude or p hase as would result from an eccentric perturber (e.g., run A13 of Figure 2 



m 



Borkovits et al.ll2003l ). To fit A ttv , we have M 2 « O.25M . The radial velocity limit (at 99.73% 



Planet b would also speed up and down twice per orbit due to a static tidal term, but this is not observable, as 
transits occur only once per orbit. 
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confidence) on circular orbits at P2 = 2Ptt v is -fT < 22 m s _1 , so the orbit would need to be inclined 
by %2 < 0.2° to the plane of the sky. Given an isotropic prior, this configuration has a probability 
of f — cos 22 = 6.0 x 10 -6 , i.e. it is too finely tuned toward face-on to be plausible. 

Having thus demonstrated that the TTV signal cannot be created by any known mechanism 
other than a second planet, we interpret the transit timing effect of planet b to be due to a second 
planet in the system, which we call Kepler-19c. 



6.1.2. Possible orbits of planet c 

In this subsection we discuss the possible orbits of planet c, consistent with the TTV data 
for planet b. We take as constraints a sinusoidal transit timing signal as well as radial velocity 
upper limit of M c < 0.44M Jup x {P c /P b ) 1/3 (99.73% confidence limit on a second orbit of arbitrary 
eccentricity). This assumes sinz c > 0.5, as we consider a mutual inclination i m > 60° to be 
physically implausible^. These planetary scenarios fall broadly into five categories: orbits with 
the period of the TTV signal, resonant perturbers, orbits near first-order mean-motion resonances, 
orbits near higher-order resonances (in which category we assign an upper limit on the mass of 
the perturber), and satellite scenarios. The latter three are favored, under circumstances that we 
elucidate below. 

Our first consideration is the possibility that planet c could have a period P c = P ttv = 316 days 
and a large eccentricity, such that it produces a time-variable tidal force on t he Kepler-19 / Kepler - 



19 b pair, accounting for the TTVs (jBorkovits et al 



2003 



Agol et al.ll2005l : iBorkovits et al. 



2011 



However, radial velocity constraints require its mass is M c < 1.4Mj up i ter , so t o generate a TTV 
signal with amplitude A ttv = 5 min, its eccentricity would need to be > 0.99 (jAgol et al.ll2005l ). 
This would be unstable with respect to planet b because the orbits would cross. Moreover, a smaller 
amplitude could be generated at 0.5 < e c < 0.9, but the signal would have a saw-toothed shape, 
which is inconsistent with the measured transit times. Therefore, any putative planet with period 
316 days that is massive enough to create the TTV signal would be inconsistent with the observed 
radial velocities. 

Another possibility is that the TTV curve of b is driven by a resonant perturb er. In this case, 
the amplitude can be substantial, even for a low-mass perturber (jAgol et al.ll2005l ). 



St 



Pb M c 
p M b + M c 



(6) 



where p refers to the period ratio of the transiting planet to the perturbing planet p/q, fi = 
max(Mb, M c )/M+, and 5t is the amplitude of the TTV signal. The libration (and TTV) period is of 



2 A large mutual inclination would likely drive large-ampli tude eccentricity cycles jKozaj 1962h in the innermost 
planet, which would in turn trigger rapid tidal orbital decay I Fabrvckv fc Tremainel 20071 : Mardline 2010h . 
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~ e 1 l 2 [i l l 2 P/p ( Agol et alj|2005l ). This consideration determines whether such a reso- 



order P\^ 

nance is a possible site for planet c, as the mass of c must be large enough so that this libration period 
equals the rather short TTV period. For low eccentricities (eh, e c < 0.1), the mass of the perturbing 
planet must be ~ 1 Jupiter mass, which is ruled out by the radial velocities. For high eccentricities, 
lower masses are allowed, but then the region of stability becomes more constrained, such that 
the proposed system would need to be finely tuned. Orbits near the Lagrange points of b (the 
1:1 resonance: Ford fc Holmanl 120071 ) fail for the same reason: P hh ~ P b a/27/4 x M+/(M b + M c ) 
(jFord fc Holmanll2007l ). which means M c ~ 3Mj up , in violation of the radial velocity constraint. 
For these reasons, we do not favor resonant orbits for planet c. 

Next, we consider perturbers near period commensurabilities. If two planet s are near to, but 



offset from, a period cq mmensurability, they can generate a large TTV signal ([Agol et al 



2005 



Hol man Murravi 120051 ) . as in the interaction between planets Kepler- 11 b and c (jLissauer et al 
201 ll ). In this scenario, the orbital frequency ratio of the known planet to the perturbing planet 



Pb/P c , is slightly offset from a ratio of integers, p/q. Then, as long as the planets' eccentricities are 
moderate, the transit timing signature will have a dominant frequency equal to: 



yPttv = \P/Pb - q/Pc 



(7) 



The absolute value sign above means we can postulate a perturber on one side or the other of each 
resonance, that is responsible for the TTV signal. The strength of the resonance must be finely 
tuned, so that a 5 minute amplitude signal is possible, yet the radial velocity amplitude also lies 
within the observed upper limit. If the resonance is spaced too closel y to the transiting planet, the 



planets will generat e timing variations on the conjunction timescale (jNesvornv &: Morbidelli 



2008 



Holman et al.ll2Q10f) which are not observed, or else the two planets will not be stable with respect 
to one another (|Wisdomlll980l ). 



The absolute value sign above means we can postulate a perturber on one side or the other 
of each resonance, that is responsible for the TTV signal. The resonance must be strong enough 
that a perturber small enough to fulfill the radial velocity constraints can produce the 5 minute 
amplitude timing signal. However, if the resonance is spaced too closely t o the transiting planet, 
the planets will generate timing variations on the conjunction timescale (jNesvornv fc Morbidelli 
2008j; iHolman et alJl2010 | ) which are n ot observed, or else the two planets will not be stable with 
respect to one another (IWisdoml ll980). 



For first-order resonances in which q = p± 1. 1 Agol et al.l (|2005l ) developed an order of magnitude 
estimate for TTV signals which depends on the fractional offset from the resonance e = |1 — 
qPb/(pPc)\ and the mass ratio jj, = max(M^, M c )/M*: 



6t h 



{M c /M*)e- x P b , e > ^ 2 

min(M c /M fe , l)// 2 e - 3 A, e < fi 1 / 2 . 



(8) 



In the current case, presuming planetary masses of < 10M®, the upper expression of Equation [S] 
holds for p < 5. In particular, for exterior first-order resonances q:p = (2:1, 3:2, 4:3), to produce 
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-Attv = 5 min we have M c ~ (4, 2, 1)M®. Of course, these values are only good to order-of- 
magnitude, but they demonstrate that a reasonable planetary mass just offset from first-order 
resonances can indeed cause the observed TTV signal. 

A planet near a second-order (e.g., 3:1, 5:3) or higher-order (e.g., 4:1, 5:1) resonance can also be 
responsible for the TTV signal. However, in these cases the strengths of the resonances are smaller, 
and they depend on a higher power of the eccentricities. Therefore, planetary-mass perturbers that 
satisfy the radial velocity constraints might need to have substantial eccentricities to match Attv 
For n : 1 resonances of exterior perturbers with n 3> 1, e c must be large, and the timing signal 
would be a series of constant-period segments with kicks at the outer body's periastron passage. 
The scaling relation of the radial velocity limit quoted above breaks down at high eccentricity, as 
periaston passages can appear as spikes which fall in data gaps, but we still wish to limit M c . 
We therefore set 99.73% confidence limits on M c sini c by (1) introducing a 4.8 M®, circular and 
edge-on planet b, and (2) sampling the orbit of planet c on a grid with P c drawn from equation [7] 
for all n : 1 resonances with n = 3 to 16 (28 cases); e c drawn from 0.1, 0.2 ... until the orbit 
crossed with planet b's; u drawn from 0°, 45°, ... 315°; To iC drawn from 10 values uniformly spaced 
between BJD 2455200 and 2455200 + P C ; and M c spaced logarithmically by 0.25 dex from 0.32M Jup 
to 18Mj up . The total grid sampled 150120 trials, and not one of the cases at or above 5.6Mj up fit 
the radial velocities (allowing for a constant offset). Apparently, although a high- mass perturber 
can be fine-tuned to not induce large radial velocity at the times of the data, its mass is limited 
even at arbitrary eccentricity. For concreteness, let us describe the end-member of this set of n:l 
resonances. It has P c = 153.2 days ~ P^ v /2, and the O — C signal would be a zig-zag, only 
marginally consistent with the apparently sinusoidal shape. At each periastron passage of planet c, 
Ph would need to change by 0.6 minutes to match A ttv - A consistent set of parameters according to 
(jHolman &; Murravll2005l . eqn. 2) is 1 Mj up it e r and e c > 0.8. The periastron of planet c would thus 
be at < 1.3a&, only marginally stable with respect to the inner planet. This scenario is also barely 
allowable according to the radial velocities. Therefore we set an upper limit of P c < 160 days. 
Again, this corresponds to M c < 6Mj up i ter , according to our grid search. 

The final possibility that we consider is a satellite orbit, which could also cause TTVs (e.g., 
Kippinell2009l ). The amplitude of 5 minutes translates to a displacement along the orbit of 2.1R p . 
A prograde orbit for such a satellite would lie be tween R p and ~ 7R P (0.4 of the Hill sphere, if 
Mb ~ 5MEarth) to be stable (jDomingos et al.ll2006l ). Since the putative satellite lies so close to the 
planet in this case, its mass would have to be > 0.6Mb cause the TTVs observed for the transiting 
planet. Therefore, it would probably be big enough to be seen in transit. We examined each 
transit by eye, to see if any deviated significantly from the single-planet mod e l, as mutual events 
of the co-orbiting planets would cause shallower transits (jSzabo et al.l 120061 : ISimon et al.l 120071 : 
Ragozzine fc Holmanll2010i ). but we found no features of interest. Furthermore, in this scenario, 
the b-c mutual orbital period would need to be near-resonant with the pair's orbital period around 
the star, so that the TTV signal aliases to the long P ttv = 316 day signal. We find this scenario 
unlikely. 



-37- 



A retro grade orbit fo r a satellite could be sta ble to much larger distances, even beyond the 
Hill sphere (|Jacksonlll913l ; IShen &: Tremaind 120081 ) . In fact, planet c could follow an independent 
Keplerian orbit which resona tes with planet b, keeping their periods and orbital phases the same 
(jLaughlin fc Chambe rs 2002). The allowable M c could be much lower in this case, down to of order 
O.lMEarth if the eccentricity is e c ~ 0.1. The TTV could be caused by either the resonant libration 
or by the two orbits precessing together at a swift rate. If the latter, then as in Section I6.1.1[ e& 
would be ~ 0.0023; unlike in that section, an apsidal motion period of 316 days is plausible because 
of the proximity of the perturber. 

We compiled examples of our favored orbits for planet c into Figure [331 Foremost are orbits 
nearby first-order mean motion resonances, which can fit the TTVs with masses which comfortably 
obey the RV constraints, and have near-circular orbits. Second, we consider orbits nearby higher 
order resonances to be possible, particularly if eccentricities are non-zero. Third, a retrograde 
satellite is a possibility, but we recognize this option as an exotic one. There is clearly profound 
degeneracy of interpretation. Such degeneracies are intrinsic to the T TV method of planet discovery , 
in the case that the signal is well-characterized by a single sinusoid ([N esvorny & Morbidclli 2008) 
and the radial velocity data cannot pinpoint the perturber ( Meschiari Sz Laughlin 201ol ). 



6.2. Constraints on Transits of Perturber 

The Kepler team searche s for transit signatures using the Transiting Planet Search pipeline 



module (| Jenkins et al.l l2010bl ). This has already been applied to the Kepler-19 light curve, and 
it identified transits of Kepler-19b. We analyzed the size of the planet we could have detected in 
the Kepler photometry by injecting signals at random phases and at varying planetary radii and 
attempting to blindly recover them. While the mass of Kepler-19c is highly uncertain, we can 
rule out the transits of at least some putative perturbers, as we describe here. The degeneracy 
of interpretation of the TTV signal means that we cannot address an exhaustive list of potential 
planets, but we comment here on the representative cases shown in Figure [Til The perturbing 
planet Kepler-19c may also not transit, in which case its orbit may be significantly misaligned from 
the transiting planet. Because Kepler-19b resides in a near-equatorial orbit (with i >88.62 with 
3a confidence), the orbit of Kepler-19c may also be misaligned by at least 0.27°, or 3.6a in the 
exterior 2:1 resonance, for example, if it does not transit. In contrast, in the interior 1:2 resonance 
at 4.57 days, the planet would have to have i <85.7 to avoid transit. The planets would have to be 
misaligned by nearly 3 degrees in that case, if the orbits are circular. 

First, we consider the orbits of 19c in second-order mean-motion or higher-order resonances. 
The planets depicted in Figure [T4l range from 1.6 to 13.8 Mg. We evaluate a minimum physical ra- 
dius for the 1.6 planet, assu ming maximum colli sional stripping of the planet during formation, 



from the relationship derived bv lMarcus et al.l (|2010i ). At the maximum possible iron fraction, a 1.6 



planet would have a radius of 1 R®. At a period of 30 days (which range encloses mean-motion 
resonances up to 3:1 with Kepler-19c) we achieve 95% completeness at 1.0 R®. At 40 days, which 
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Fig. 14. — Possible orbits for Kepler-19c. Orbits near first-order mean-motion resonances (eq. [7]) 
may fit the TTV signal and RV constraints even on circular orbits; shown are P c = 6.129 days 
(bright green) and 6.256 days (blue), flanking the interior 2:3 resonance, and P c = 18.033 days 
(olive green), next to the exterior 2:1 resonance. Orbits near higher-order resonances likely need 
eccentricity to produce the TTV signal with masses low enough to satisfy the radial velocities; 
shown are P c = 3.065 days (pink) near the interior 1:3 resonance, P c = 15.326 days (purple) near 
the exterior 5:3 resonance, P c = 27.036 days (aqua blue) near the exterior 3:1 resonance, and 
P c = 38.310 days (brown) near the exterior 4:1 resonance. Finally, a co-orbital planet (or distant 
retrograde satellite) is shown at P c = 9.287 days (red) - in such an orbit very small masses are 
possible, so the dot representing the planet is drawn small. Other possible orbits are within mean 
motion resonances (including the 1:1 at Lagrange points), or a prograde satellite, but these are 
disfavored (see text). 



includes the example perturber in Figure [T4"l in the 4:1 mean- motion resonance, we would still have 
detected a 1.0 R® planet at 90% of phases. 

If the planet Kepler-19c were instead coorbital with 19b, or if it resides in a satellite orbit, its 
mass could be much smaller, as described in the previous section: this mass could be as smaller 
as 0.1 M ffi ( e quiva lent to the mass of Mars) if its eccentricity were equal to 0.1. The models of 
Seager et al.l (|2007l ) show that a 0.1 Mq, planet could be as small as 0.4 R@ if it comprised 70% 
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iro n and 30% silicate b y mass (which is plausible, given the maximum iron fractions determined 
by iMarcus et al.l l2010l ). Its predicted transit depth would be 20 ppm, which might be barely 
detectable when compared to the error bar on the transit depth of Kepler-19b (with a necessarily 
similar orbital period) of 6 ppm. However, even if such a planet transited, its detection would be 
extremely challenging. 

However, though the representative cases (with the exception of the coorbital scenario) depicted 
in Figure [H] would all have been readily detectable, we note that the perturbing planet could easily 
be smaller than 1 in the mean- motion or higher-order resonance cases, in which case it might 
have eluded detection even in orbits with periods shorter than 30 days. However, we note that we 
achieve 95% completeness for 0.7 -R^-sized planets up to 10 days, so such a world would have to be 
less massi ve than 0.7 M m for P < 10 days, referring again to the maximum iron fraction models 
derived by Marcus et al.l (|2010l ) . 



6.3. Search for Secondary Eclipse of Kepler-19b 

If we assume the planet reradiates isotropically the energy it receives from its star, then the 
equilibrium temperature of Kepler-19b is given by: 

T p = (l-A B )V 4 n^, (9) 

where T* is the temperature of the star, a/R* is the orbital radius to stellar radius ratio, and Ab 
is the albedo of the planet. If we assume a Bond albedo Ab of 0.3 and ignore atmospheric effects 
in order to obtain a rough estimate for the equilibrium temperature, we may employ the MCMC 
chain of a/R± (and the corresponding values for T e g of the star identified from the nearest stellar 
isochrone, per the analysis in Section 3.2) to find the allowable range of planetary temperature. 
The range that encompasses 68% of realizations of planetary temperature is Tp = 770 ± 10 K, 
rounded to the nearest 10 Kelvin. 

There are two contributing sources to the occultation depth, namely the reflected starlight, 
and the emitted light from the planet itself. Both of these are dependent upon the unknown albedo 
of the planet. Assuming again the Bond albedo of 0.3, the expected depth due to reflected light is 
given by 5 re f = AsiRp/a) 2 ; this value lies between 3 and 4 parts in 10 7 . The expected occultation 
depth due to the emitted light of the planet is given by 5 em = (Rp/R*) 2 ■ B\{T p ) / B\(T+) . Using a 
wavelength of 700 nm (in the middle of the Kepler bandpass) to estimate B\(T), S em is of order 
10~ 13 and so contributes negligibly to the expected eclipse depth. 

To assess our sensitivity to the secondary eclipse, we fit a line to either side of each expected 
eclipse, and divided this line from the data during eclipse (centered on a phase of 0.5 and spanning 
2.5 eclipse durations on either side), in a manner similar to the method we employed for data in 
transit. We then evaluated a model with epoch and duration set by the transit parameters given in 
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Table [U but with variable depth, from to 20 ppm (considerably larger than the expected eclipse 
depth), and compared to the phase- folded light curve. We neglect the possibility that the transiting 
planet resides in an eccentric orbit, and thus that the secondary eclipse may not occur at a phase of 
0.5. We find that our sensitivity is yet too low to detect the secondary eclipse at any phase. Values 
as high as 9 ppm are statistically indistinguishable from a depth of zero (furnishing a \ 2 difference 
less than 9). A planet with a maximal Bond albedo of 1 would produce a decrement of 1.2 ppm, 
which is still considerably below detectability with Kepler. Figure [15] shows the x 2 improvement 
associated with adding an eclipse of variable depth at a phase of 0.5. The dotted line shows the 
expected depth for an extreme albedo of 1, which is indistinguishable from a flat line. 

1 .00006 F ' ' ' ' ' r ~f "~T 

1.00004 '- t t T 

1 .00002 '- T T T T T | T T t -. 

2 1.00000 i-L^4*TT""-4^^i-4tL-— — ^■■^"■■^™~r 

U. 1 I 1 ^ I T 1 T t 

0.99998 ]■ T 1 1 | 1 1 1 } | 1 t 1 1 1 1 T 1 \ 

0.99996 r 1 [ 1 

0.99994 F 1 : 

-3-2-10 1 2 3 

Phase [hours] 




2.0X10" 6 4.0X10" 6 6.0X10" 6 8.0x10" 6 1.0x10" 5 1.2x10" 5 

Depth of Secondary 

Fig. 15. — The Kepler-19 light curve at a phase of 0.5. The best eclipse light curve is shown 
overplotted, but this solution is statistically indistinguishable from a flat line. The expected depth 
for planet with an albedo of one is shown by the dashed line. 



6.4. Composition of Kepler-19b 



While we cannot estimate the mean density of Kepler-19b without a measurement of its mass, 
we can still dra w meaningful concl usions about its composition from the upper limit mass value, 
as was done bv lFressin et al.l ()201ll ) in the similar case of Kepler-lOc. We first address whether we 



-41 - 



can rule out solid compositions at the highest density: a planet made of pure iron at a radius of 
2.2 Rq would have a mass of 100 M® (jSeager et al.ll2007l ). However, such a high fraction of iron is 
unphysical, even with maximal collisional stripping during the planet's formation; a 2.2 R& , planet 
with the largest possible iron fraction would have a mass of 30 M® (jMarcus et al.ll2010i ). This 
maximum density is ruled out with 95% confidence by the mass upper limit from radial velocities 
of 20.3 M®. A planet composed o f pure silicate at the measured radius, however, would have a 
mass of 15 M® (jSeager et al.1 120071 ). which lies within the allowable mass range for Kepler-19b. In 
cont rast, a 2.2 i?® plan et with a homogeneous composition of water ice would have a mass of 4.5 



M® (jSeager et al.l 120071 ): mixtures of water ice and silicate in any fraction are therefore consistent 



with the measured mass upper limit. We consider also whether a substantial H/H e envelope is 



possib le for Kepler- 19b. This scenario brackets the lower range of possible densities. iRogers et al 



(|201ll ) present theoretical models for planets in the radius range of 2-6 it!® and the temperature 
range 500-1000 K (in which sample Kepler-19b, with radius of 2.2 i?® and temperature near 700 K, 
is included), given different formation histories. If Kepler-19b formed by a nucleated core-accretion 
scenario beyond the snow line, a core of ice and rock surrounded by an H/He envelope would not be 
tenable even at the cooler equilibrium temperature of 500 K: such an envelope would have been lost 
in a timescale of <1 Gyr. An outgassed hydrogen envelope, by comparison, is a plausible scenario 
over a timescale greater than 1 G yr, although the m ass fraction of such an envelope would be less 



than 0.01 the mass of the planet ([Rogers et al.l 12011 



6.5. Future Prospects 

In light of the dynamical study of Section 6.1, which combined transit time variations and radial 
velocities to characterize a perturber in the Kepler-19 system, we may ask whether we can expect 
to measure the precise orbital parameters of planet 19c. The two main issues causing ambiguity 
among the perturber scenarios are: (i) the transit variations are consistent with a smooth sinusoid 
down to the noise l evel, with no additional hints of the perturber's identity such as "chopping" 



(jHolman et alJl2010l ) on the conjunction timescale, (ii) the radial velocities, while essential to ruling 
out massive, non-planetary perturbers, simply do not have the precision sufficient to distinguish 
between the various planetary scenarios. The former issue might be resolved with much more data. 
The Kepler collaboration intends to keep Kepler-19 on its short cadence mode (demonstrated here 
to result in superior timing accuracy) for the remainder of the mission. Furthermore, in numerical 
simulations, chopping signals can wax and wane over secular cycles, so a detection of this effect is 
possible in future quarters. The latter issue requires more radial velocity data to resolve, which 
could allow us to further address the degeneracy of interpretation of the TTV signal, or could point 
to even more planets that are not clearly detected by transit timing variations. 
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with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California 
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Table 2. Transit Times for Kepler-19b From Q0-Q8 



Transit Number Transit Time Predicted Time 0-C — la +la 

[BJD-2450000] (from linear ephemeris) [min] [min] [min] 



1 


4959.70744 


4959.70605 


2.0 


2.6 


2.7 


2 


4968.98895 


4968.99305 


-5.9 


3.4 


5.6 


3 


4978.27949 


4978.28004 


-0.8 


3.6 


3.2 


4 


4987.56280 


4987.56704 


-6.1 


3.3 


3.2 


5 


4996.85396 


4996.85403 


-0.10 


3.8 


3.5 


6 


5006.13936 


5006.14102 


-2.4 


3.4 


4.2 


8 


5024.71126 


5024.71501 


-5.4 


3.2 


3.3 


9 


5033.99971 


5034.00201 


-3.3 


4.0 


3.9 


10 


5043.28879 


5043.28900 


-0.30 


3.2 


3.5 


11 


5052.57509 


5052.57599 


-1.3 


3.9 


4.0 


12 


5061.85993 


5061.86299 


-4.4 


3.2 


3.3 


13 


5071.14832 


5071.14998 


-2.4 


3.5 


3.4 


14 


5080.43524 


5080.43698 


-2.5 


3.3 


4.5 


15 


5089.72744 


5089.72397 


5.0 


3.2 


2.9 


16 


5099.00860 


5099.01096 


-3.4 


2.2 


1.7 


17 


5108.29463 


5108.29796 


-4.8 


1.7 


1.5 


18 


5117.58342 


5117.58495 


-2.2 


2.0 


1.7 


19 


5126.87077 


5126.87195 


-1.7 


2.8 


2.0 


20 


5136.15957 


5136.15894 


0.9 


4.1 


1.5 


21 


5145.44635 


5145.44593 


0.6 


1.8 


2.0 


23 


5164.02277 


5164.01992 


4.1 


2.0 


0.9 


24 


5173.30976 


5173.30692 


4.1 


1.5 


1.3 


26 


5201.17269 


5201.16790 


6.9 


1.5 


1.3 


27 


5210.45913 


5210.45489 


6.1 


3.9 


3.1 


28 


5219.74591 


5219.74189 


5.8 


1.4 


1.4 


29 


5229.03319 


5229.02888 


6.2 


1.6 


2.3 


30 


5238.31921 


5238.31587 


4.8 


1.3 


1.3 


31 


5247.60551 


5247.60287 


3.8 


1.4 


1.4 


32 


5256.89243 


5256.88986 


3.7 


1.5 


1.8 


33 


5266.18033 


5266.17686 


5.0 


1.5 


2.0 


34 


5284.75126 


5284.75084 


0.6 


2.6 


1.6 



Table 2 — Continued 



Transit Number Transit Time Predicted Time O-C — la +la 

[BJD-2450000] (from linear ephemeris) [min] [min] [min] 



3o 


5294.03582 


ron/i no to \ 

5294.03784 


-2.9 


1.0 


1.9 


37 


5303.32205 


5303.32483 


-4.0 


1.7 


2.2 


38 


p* o -i o r> -\ r\-\ r> 

5312.61016 


5312.61183 


-2.4 


1.6 


1.7 


on 

39 


5321.89569 


p* o o 1 onoon 

5321.89882 


A P* 

-4.5 


2.0 


1.3 


40 


5331.18179 


p* n o ~i i o r* o 1 

5331.18581 


p* o 

-5.8 


1.5 


1.3 


A 1 

41 


ro /in .4 ^ o p* T 

5340.46857 


ro /in t Tooi 

5340.47281 


-6.1 


2.6 


1 T 

1.7 


A ri 

42 


rn «n TP* P* TT 

5349.75577 


5349.75980 


-5.8 


O T 

8.7 


1.6 


A O 

43 


5359.04291 


rnrn a «/>on 

5359.04680 


-5.6 


1.6 


1 T 

1.7 


A A 

44 


5368.33171 


p* o r* o oooTn 

5368.33379 


-3.0 


1.5 


2.3 


45 


p*ott / ' i <"i • > i"i 

5377.61939 


roTT o o t o 

5377.62078 


-2.0 


1.8 


1.5 


4b 


5386.90368 


5386.90778 


-5.9 


2.0 


1.8 


A >"7 

47 


5396.19179 


5396.19477 


A O 

-4.3 


C\ A 


1 T 

1.7 


48 


5405.47968 


5405.48177 


-3.0 


2.4 


1.5 


49 


5414.76688 


5414.76876 


-2.7 


2.2 


1.5 


50 


5424.05381 


5424.05575 


-2.8 


3.2 


2.4 


51 


p* J o o o ,1 o T r* 

5433.34275 


5433.34275 


0.0 


1.8 


2.1 


52 


5442.62821 


5442.62974 


-2.2 


3.2 


3.2 


53 


5451.91715 


5451.91674 


0.6 


2.0 


1.4 


54 


5461.20241 


5461.20373 


-1.9 


1.7 


3.8 


55 


5470.49253 


5470.49072 


2.6 


1.8 


1.9 


56 


5479.77723 


5479.77772 


-0.7 


1.3 


1.4 


57 


5489.06589 


5489.06471 


1.7 


2.4 


3.1 


58 


5498.35615 


5498.35171 


6.4 


1.8 


2.1 


59 


5507.64127 


5507.63870 


3.7 


1.3 


1.2 


60 


5516.92944 


5516.92569 


5.4 


1.5 


1.8 


61 


5526.21866 


5526.21269 


8.6 


2.1 


1.8 


62 


5535.49975 


5535.49968 


0.1 


1.4 


1.7 


63 


5544.78855 


5544.78668 


2.7 


1.3 


1.2 


64 


5572.64863 


5572.64766 


1.4 


1.7 


1.3 


67 


5581.93437 


5581.93465 


-0.4 


1.4 


1.8 


68 


5591.22095 


5591.22165 


-1.0 


1.2 


1.2 



-50- 



Table 2 — Continued 



Transit Number 


Transit Time 


Predicted Time 


O-C 


-1(7 


+ 1(7 




[BJD-2450000] 


(from linear ephemeris) 


[min] 


[min] 


[min] 


69 


5600.50788 


5600.50864 


-1.1 


1.3 


1.4 


70 


5609.79355 


5609.79563 


-3.0 


1.3 


1.6 


71 


5619.07832 


5619.08263 


-6.2 


1.4 


1.7 


72 


5628.36705 


5628.36962 


-3.7 


1.2 


1.2 



Note. — Transits with numbers < 16 were gathered at long cadence (with an exposure 
time of 29.5 minutes), while transits with numbers > 16 were gathered at short cadence 
(with an exposure time of 58.8 seconds). 



